Astrophysics and Space Science 
DOI 10.1007/s»«»«-»»«-»«»«-» 
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Abstract An overview of selected topical problems on 
modelling oscillation properties in solar-like stars is pre- 
sented. High-quality oscillation data from both space- 
borne intensity observations and ground-based spectro- 
scopic measurements provide first tests of the still-ill- 
understood, superficial layers in distant stars. Empha- 
sis will be given to modelling the pulsation dynamics 
of the stellar surface layers, the stochastic excitation 
processes and the associated dynamics of the turbulent 
fluxes of heat and momentum. 

Keywords asteroseismology; convection; pulsation 
mode physics; stellar structure. 



1 Introduction 



With the high-qua lity photometric data from the Ke 



pier satellite (e.g. IChristensen-Dalsgaard et al.1 I2007T ) 



and spectroscopic data from ground-based observation 
camp aigns, and later also fro m the Danish SONG net- 



work (jGrundahl et al] 120071 ). we shall be able to ad- 



dress more carefully many of the current problems 
in modelling pulsation properties in solar- like stars. 
In solar-like stars the p-mode lifetimes and ampli- 
tudes are crucially affected by the processes that take 
place in the outer convectively unstable stellar layers. 
A proper modelling of the dynamics of the convec- 
tive heat and momentum transfer is therefore essen- 
tial. Recent 3D numerical si mulations of the largest 
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2008) have 


been proven to be extremely useful for 


calibrating 
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and testing semi-analytical models for convection and 
stochastic excitation. However, the high-Reynolds- 
number (and low-Prandtl-number) turbulent convec- 
tion in stars still prohibits today's simulations from 
resolving all the required scales of stellar turbulence. 
Consequently such simulations have to use sub-grid- 
scale models, which may lead to different results (e.g. 



Jacoutot et al.l 120081 ) . Because this situation will not 
change in the near future, we still need analytical mod- 
els for describing the convection and pulsation dynam- 
ics in stars. 

Sections 2 and 3 will discuss selected problems of 
our current understanding of nonadiabatic pulsation 
dynamics in the Sun and in the hotter F5 star Pro- 
cyon. Although we can reasonably well reproduce the 
observed pulsation prope rties in the Sun, th e recent 
Procyon observations (e.g. Arentoft et al-lEooil and ref- 
erences therein) have revealed serious problems in our 
models for estimating the oscillation amplitudes in stars 
hotter than the Sun. In Section 4, therefore, I shall 
address the problem of selecting a proper temporal tur- 
bulence spectrum for modelling the stochastic energy 
supply rate for acoustic modes. 



2 Mode parameters 

In solar-like stars all possible modes of oscillation are 
stable; thus, if a given oscillation mode is somehow ex- 
cited, its amplitude will decay over a finite time, typi- 
cally of the order of days to months, the inverse of which 
is the damping rate r\. The oscillation power spectrum 
can be described in terms of an ensemble of intrinsically 
damped, stochastically driven, simple-harmonic oscilla- 
tors, provided that the background equilibrium state of 
the star were independent of time. In that case the 
mode profile is essentially Lorentzian, and the intrin- 
sic damping rates of the modes could then be deter- 
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Fig. 1 Top: Calculated amplitude ratios (see equation 
Q) as a function of height in a solar model for modes with 
different frequency values. Bottom: Theoretical amplitude 
ratios (surface luminosity perturbation over velocity) for a 
solar model compared with observations by Sc hriiver et al.l 
(jl99ll ). Computed results are depicted at different heights 
above the photosphere (/i=0km at T = T c g). The thick, 
solid curve indicates a ru nning-mean average of the data 
(from iHoudek et aDll999l ) 



mined observationally from measurements of the pulsa- 
tion linewidths. The other fundamental property of the 
observed oscillation power spectrum is the height, H, 
of a single peak in the Fourier spectrum. The observed 
velocity signal v(t) — d£/di (where £(i) is the surface 
displacement of the damped, stochastically driven, har- 
monic oscillator, and t is time) can then be related to 
the mode height H by taking the Fourier transforma- 
tion of the harmonic oscillator followed by an integra- 
tion over frequency to obtain the total mean energy 
£ in a particular pulsation mode with inertia / (e.g. 



Chaplin et al. 200 



pulsation 
3 iHoudekl 



20061 ). For Ti] > 1, where 



T is observing time, the squared surface rms velocity is 
then given by 
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Fig. 2 Calculated amplitude ratios (see Eq. Q) for a 
model of Procyon A are compared with observations by Ar- 
entoft (2008; horizontal dotted line). Theoretical results are 
shown for a scaled VAL-C atmosphere (black, thin curves) 
and for an Eddington atmosphere (red, thick curves). Top: 
The theoretical amplitude ratios are shown as a function of 
height and for three different pulsation modes. The frequen- 
cies of the three pulsation modes are indicated. Bottom: 
The theoretical amplitude ratios are shown as a function 
of frequency at three different heights in the stellar atmo- 
sphere. The heights above the photosphere h = km are 
indicated 

where P is the energy supply rate in ergs -1 , and H 
is given in units of cm 2 s _2 Hz _1 . The height H is the 
maximum of the discrete power, and is obtained from 
integrating the power spectral density over a frequency 
bin: 



H= / \V{v)\ 2 du. 

Ju-S/2 



(2) 



where V(y) is the Fourier transform of v(t), v is cyclic 
frequency, and 6 = 1/2T, is the frequency bin deter- 
mined by the observation time T. It is therefore not the 
total integrated power, V 2 , that is obs erved directly. 



but r ather the power spectral density (jChaplin et al 
20051) . 
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3 Pulsation amplitude ratios 

Linearized pulsation equations for nonadiabatic radia l 
oscillations can be presented as (e.g. iBalmfort h 1992a): 



_d_ /Sp 
dm \ p 
_d_ ?Sr 
dm \ r 
_d_ (5L\ 
dm \ L J 



(5r ST 6p Spt <5$ 

1 f 3 $r Sp_ 

Awr 3 p \ r p 

c p T f ST 6p 



(3) 



where 5 is the Lagrangian perturbation operator, and 
for simplicity the right hand side of the perturbed mo- 
mentum equation is formally expressed by the func- 
tion /. Equations ([3|), together with an equation relat- 
ing the heat flux to the stratification of the star (the full 
set of equations can be found in, e.g., Balmforth 1992a), 
are solved subject to boundary conditions to obtain the 
eigenfunctions and the complex angular eigenfrequency 
LOi = uj^+irji of the mode i, where tu T i is the (real) pulsa- 
tion frequency and r\i the linear damping rate in (s -1 ). 
The turbulent flux perturbations of heat and momen- 
tum, SL C and Sp t , and the fluctuating anisotropy fac- 
tor S& are obtained from t he nonlocal, tim e-dependent 
convection formulation bv lGoug hi (Il977al lbh. 

From the linearized nonadiabatic pulsation equa- 
tions ([3]) theoretical intensity- velocity amplitude ratios 



SL/L 



(4) 



AV uj^r Sr/r 

can be compared with observations, without the need 
of a specific excitation model and all its uncertainties 
in describing the turbulence spectrum. 

In the top panel of Fig. [T] the theoretical amplitude 
ratios (equation Q) of a solar model are plotted as a 
function of height for several radial pulsation modes. 
The square root of the mode kinetic energy per unit in- 
crement of r, which is proportional to rp 1 ' 2 Sr, increases 
rather slowly with height; the density p, however, de- 
creases very rapidly and consequently the displacement 
eigenf unction Sr increases with height. This leads to 
the results shown in the upper panel of Fig. [T] The 
decrease in the amplitude ratios with height is partic- 
ularly pronounced for high-order modes for which the 
eigenfunctions vary rapidly in the evanescent outer lay- 
ers of the atmosphere. It is for that reason why solar 
velocity amplitudes from e.g., the GOLF instrument 
have larger values than the meas urements from the B i- 
SON instrument (by about 25%. iKieldsen et al.ll2005l ). 
The lower panel of Fig. Q] compares the estimated solar 
amplitude ratios (curves) with observed ratios (sym- 
bols) as a function of frequency. The model results are 



depicted for velocity amplitudes computed at different 
atmospheric levels. The observations are obtained from 
accurate irradiance measurements from the IPHIR in- 
strument of the PHOBOS 2 spacecraft with contempo- 
raneous low-degree v elocity data from the BiSON in- 
strument at Tenerife (jSchriiver et al.lll99ll ). The thick 
solid curve represents a running-mean average, with a 
width of 300 /iHz, of the observational data. The theo- 
retical ratios for h = 200 km (dashed curve) show rea- 
sonable agreement with the observations. 

In Fig.[2]model results for the F5 star Procyon A are 
compared with observ ations (horizontal dotted line) by 



Arentoft et al. ( 2008| ). Theoretical re sults are shown 



for tw o stellar atmospheres: a VAL-C (jVernazza et al 



19811 ) atmosphere scaled with the model's effective tem- 
perature T e ff (thin curves), and an Eddington atmo- 
sphere (thick curves). For both stellar atmospheres 
the agreement with the observations is less satisfactory 
than in the solar case, indicating that we may not repre- 
sent correctly the shape of the pulsation eigenfunctions. 
Consequently there is need for adopting more realisti- 
cally computed atmospheres in the equilibrium models, 
particularly for stars with much higher surface temper- 
atures than the Sun. It should, however, be mentioned 
that the current photometric observations are still un- 
certain. 



4 Stochastic excitation and turbulent spectra 

In solar-like stars the oscillations are driven by the vig- 
orous turbulent convection in the very outer surface 
layers. The turbulent fluid motion generates acoustic 
waves in a broad frequency range, which excite a large 
number of global p modes (of the order of 10 million 
p modes in the solar case), and poss ibly also global g 
modes (e.g. Appourchaux et a"l1l2009h . In the past sev- 
eral stochastic excitati on models were proposed (for re- 



cent reviews see, e.g. lHoudekll2006l : lAppourchaux et al 
2009h . In general all reported excitation models re- 



produce reasonably well the observed oscillation ampli- 
tudes for solar-like stars that are similar or cooler than 
the Sun. However, for stars that are somewhat hot- 
ter than the Sun, theoretical predictions overestimate 
the pulsation amplitudes by up to a factor of about 
four, such as for the F5 star Procyon A. Additional to 
the uncertainties in modelling the shape of the pulsa- 
tion eigenfunctions in the atmospheric stellar layers (see 
Secti on El), modell ing of the mode damping rates (see 



e.g., lHoudekll2006f ) and the turbulent energy spectrum 



play also a crucial role in estimating the mode ampli- 
tudes. In this paper I shall only address the problem 
of modelling the turbulent spectrum for estimating the 
energy supply rate P. 
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Fig. 3 Comparison of various analytical frequency factors 
fi(aj, Tk',r) [Eq. pil) -(!13 p ] that are used in stochastic excita- 
tion models 



Fig. 4 Effect of varying the correlation parameter A 
[Eq. (JBJ] ° n the Lorentzian frequency factor t^; r) 

[Eq. (H3j)] 



4.1 Stochastic excitation model 



A revi ew on this topic was recently given bv lAppourchaux 
( 2009f ). Therefore I shall summarize only the most im- 
portant matters. According to Eq. ([1]) we need to model 
the (linear) damping rate r\ and the energy supply rate 
P for estimating the mode height H. In this paper 
I shall address only the modelling of the energy sup- 
ply rate. One of the problems that one faces in de- 
riving a stochastic excitation model is the separation 
of the total fluid motion into (large-scale) oscillatory 
motion (oscillation modes), with displacement £(r,t), 
and small-scale turbulent motion with the convective 
velocity field u(r,t) = (u,v,w). This separatio n is su- 
perfici ally straightforward for radial p modes (jGoughl 
Il977al ): however, for nonradial modes this separation is 
a much more difficult task and can possibly not be ac- 
complished for the very-high-degr ee p modes and also 



bulent velocity field u. Here we consider only the 
Remolds stress driving term, which dominates over 



Siyji^Jentropy driving term in the Sun and in most solar 



very- high-order g modes (see e.g. lAppourchaux et al 



2009). For simplicity we shall discuss here only radial 



p modes. 

Follow ing the basic p rinciple o f Lighthill's aco ustic 
analogy ( LighthilJ 19521 see also Crightonl 1975 ) one 
typically arranges the linearized global mode variables 
on the left and the nonlinear fluctuating terms associ- 
ated with the turbulent convection on the right of the 
equations of motion, which may then be written as 



d 2 £ „ d£ 



(5) 



in which the linear spatial differential operator C 
satisfies the (unforced) homogeneous wave equation 
££j = ujf£ for the mode eigenfunctions with fre- 
quencies uji, and J- is the nonlinear stochastic forc- 
ing and damping term that depends only on the tur- 



tvpe stars jBalmforthl[i992bl Istein h Nordlundl 



2001 



Stein et al.l 12004 iBelkacem et al.1 12006E iSamadi et al . 
20071) . It should, however, be noted that there is partial 



cancellation between the fluctuating Reynolds stress 
and entropy source terms, which could lead to smaller 
oscillation amplitudes t han adopting the Reynolds 



stress contribution alone (Osaki 1990; Stein et al. 2004 



HoudeklE)0l . Because the mode amplitudes of solar- 



type oscillations are small it is additionally assumed 
that they do not interact with the nonlinear right 
hand side of the wave equation, and consequently 
one can solve Eq. (Q by a nonsingular perturbation 
method dGoldreich fc Keelevl Il977t iBalmfortbl Il992b : 
Samadi fc Goupill l200ll IChaplin et al.ll2005h . For ra- 



dial modes only the vertical component T r of the fluc- 
tuating Reynolds stress source term F is important, 



T r (yS) — d r (pww — (pww)) . 



(G) 



which depends on the (r, r) component of a two- 
point correlation function, R rr := (ww), where angu- 
lar brackets denote an ensemble average. For incom- 
pressible, homogeneous isotropic turbulence the Fourier 
transform R rr of R rr is proportional to the turbulent 
energy spectrum function E(k,u), i.e. 



R rr oc k 2 E(k, ui) 



(7) 



( Batchelor I l953l ) . where A; is a wavenumber. Following 
Steml (|l967h~ we factorize the energy spectrum function 



into E(k, ui) 
X/kuk 



E(k£/ir)Q(uj,Tk;r), where 



Tk 



(8) 
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Fig. 5 Temporal spectrum of the normalized velocity au- 
tocorrelation (|17[) for two values of the Reynolds number 
R c , computed according to Sawford's (1991) second-order 
stochastic turbulence model 
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Fig. 6 Measured temporal velocity power spectrum in a ro- 
tating shear turbulence experiment with water (solid curve). 
The dot-dashed curve is a fit of Sawford's (1991) model, 
Eg. JTTb. to the dat a with R c = (T-l/T-j) 2 (adapted from 
iMordant et~ai1l2004l ) 



is the correlation time-scale of eddies with vertical size 
£ = Tr/k and velocity u k ', the correlation factor A is of 
order unity and accounts for uncertainties in denning 
Tfc. For statistically stationary turbulence the energy 
supply rate is then given by (see e.g., Chaplin et al. 
2005 for details) 



P=^ 
91 



with 



^rpt^) 5(w i; r)dr, (9) 



K- 2 E 2 {K)tl{LU u T k -r)&K, (10) 



where pt = (pww) is the (r, r) component of the (mean) 
Reynolds stress tensor, k = k£/ir, R is surface radius, 
£ r i is the normalized radial part of and the prod- 
uct of <!> and ^ is a factor of unity accounting for the 
anisotropy of the turbulent velocity field. The spectral 
function S accounts for contributions to P from the 
small-scale turbulence. For the normalized spatial tur- 
bulence energy spectru m E(k) it has been common to 
adopt, for example, the Kolmogorovl ( 194ll ) spectrum. 
For the frequency-dependent factor Q(cj,Tk;r), how- 
ever, which is used for evaluating the self-convolution 
Vt{uoi,T k ;r) = J Q(u},Tk]r)Cl(u>i-u),Tk;r) dco in Eq. (jTOJ) , 
no satisfactory theory exists. Various functional forms 
were proposed in the past, which lead, however, to 
rather different values for the modelled oscillation am- 



plitudes H (S amadi e ^al.1 12003 : IChaplin et alj 12005 



lAppourchaux et al.l l2009h . We shall therefore discuss 
the most commonly adopted frequency-dependent fac- 
tors for stellar turbulence spectra in the next section. 



4.2 Turbulence spectra 

Differences in the (properly normalized) spatial spec- 
trum E(k) create only minor differences in the energy 
supp l y rate P (e.g..lBalmforthlll992bl ISamadi fc Goupill 



20011 IChaplin et all 120051 ). But differences in the 
frequency-dependent (temporal) spectrum n(u>,Tk; r) 
would create rather large differences in the p-mode and 



-mode amplitudes (Samadi et al.fi2003: Chaplin et al. 


Belkacem et al. 


20091 Appourchaux et al. l2009h. 



Those differences are produced by the convection- 
oscillation interactions in the deeper convectively un- 
stable stellar layers, which are off resonance, and whose 
magnitude depends crucially on the adopted frequency 
dependence of the turbulence spectrum in the high- 
frequency tail of the cascade. We consider the following 
factors: 



(i) the Exponential factor (|Stemlll9671 ) 



r2 E (w,T fc ;r) = 



V21n2 



-r k e 



-(ujT k V2\n2/2y 



(ii) the Gaussian factor (|Steinlll967l ) 

fi G (w,T fe ;r) 



(11) 



(12) 



(hi) the Lorentzian factor (jGoug h 1977a: ISamadi et al 
2003; IChaplin et aDl2005h 



n L (w,T fc ;r) 



Tk 



7rV21n2 1 + (ur fc / 'V21n2) 2 



(13) 



Fig.[3]cornpares the dimensionless quantity Q(ui, r k ; t)t^ 1 
for the three frequency factors (fTT |) -(fl"3 |) . All three fac- 
tors are normalized according to J £l(u), r k ; r) dto = 1, 
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Fig. 7 Comparison of the temporal power spectrum of the 
vertical convective velocity w in the Sun, computed with 
Stein & Nordlund's 3D hydrocode (solid curve), with ana- 
lytical frequency factors (symbols). Results are shown for 
a depth of 250 km below the s urface and for a horizon tal 
wavenumber k — 4 Mm -1 (from [Georgobiani et al.ll2004) 



Furthermore 
l 



the 
has been 



the integral being from — oo to oo. 
full-width at half-maximum, i.e. 2\/2 In 2r t 
chosen to be the same for all three factors, in order 
to have a meaningful comparison. The Gaussian fre- 
quency factor decreases more rapidly with WTfe, than 
both the Exponential and Lorentzian factors, i.e. the 
Gaussian factor produces the smallest magnitude in the 
high-frequency tail of the spectrum. Convective eddies 
situated in the deeper convectively unstable stellar lay- 
ers have longer characteristic time scales r^, i.e. Tfc [see 
Eq. ©] i ncreases with depth in the star (see for example 



Fig. 5 of IChaplin et al.l 120051 ). The Lorentzian factor 



decreases more slowly with depth at constant frequency 
and consequently a larger fraction of the integrands of 
Eqs (J9j> and (fTOf arises from large off-resonant eddies 
situated deep in the star, whereas the Gaussian fre- 
quency factor gives less weight to the large off-resonance 
eddies. As a result of this the modelled oscillation am- 
plitudes are larger with a Lorentzian time-correlation 
function and therefore in better agreement with the 
obser vations than with a Gaussian (|Samadi et al.ll2003 . 
20071) . However, IChaplin et all ((2005J) reported that 
the Lorentzian time-correlation function leads to over- 
estimated heights H at low f requencies for solar 



modes. Also for solar g modes iBelkacem et al.l (|200 



' P 
)9) 



concluded that the best fit for reproducing the observa- 
tions was obtained with a Gaussian factor at the lowest 
and a Lorentzian at th e highest frequencies. Moreover, 
Belkacem et al.l (|2009j ) reported that the correlation pa- 
rameter A [see Eq. (jHJ] has to increase with stellar depth 
in order to reproduc e the 3D numerical simulations by 
Miesch et al.l d2008l). a result that is consistent with 
the findings bv IChaplin et alJ (|2005h . The effect of A 




10 15 
k (1 /Mm) 

Fig. 8 Exponent j3 and width A" 1 for the empirical ex- 
pression (|18p for different (mean) horizontal wavenumbers 
k and solar depths. f3 & A are obtained from fitting to the 
simulated temporal powe r spectra of the vertica l velocity 
the expression (|18p ( from [Georgobiani et alj|2004l ) 



on the time-correlation functions is demonstrated in 
Fig.@ for a Lorentzian factor. Increasing A leads to 
a more rapidly declining tail and consequently to a 
smaller contribution to the energy supply rate P from 
the off-resonant eddies situated in the deep convection 
zone. Finding the correct time-correlation function for 
stellar turbulence remains an open issue. 

From the theoretical point of view a Lorentzian fac- 
tor is a result predicted for the largest, most-energetic 
eddies by the time-depe ndent mixing-length formu- 
latio n of Goughl (1977a). Moreover, in ag reement 
with Kolmogorovl ( 1941 ) theory Sawfordl ( 199ll ) demon- 
strated that in the limit of high-Reynolds-number tur- 
bulence the two-point velocity autocorrelation, normal- 
ized by the velocity variance a 2 :— (ww) — (w)(w), is 



R(t) = (w(t)w(t + t)) /a 2 = c-l r l/ TL 
where 



TV 



R(t) dr 



(14) 



(15) 



is the Lagrangian integral time scale (see also lLegg fc Raupach 
19821) . Consequently the normalized velocity autocor- 
relation R is proportional to an exponential decrease 
with separation time r in the inertial range, and the 
velocity spectrum becomes with Eq. (|14[) 



R(t) cos(wt) dr 



(16) 
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i.e. it is of Lorentzian shape with F(u) cx u>~ 2 for li — > 
oo and at infinite Reynolds number R Q . Sawford ex- 
tended this analysis for arbitrary values of the Reynolds 
number R c := (<eA»;) 2 , where and t v are the time 
scales of the energy-containing e ddies and the Kol- 
mogorov timescale respectively ([Tennekes &: Lum lev 



1972). For this case the velocity spectrum is (jSawforc 



1991) 



F(w) 



(1 + R 



1/2- 



(R, 



1/2 



c 2 t 2 ) 2 



(1 



Rl /2 



(17) 



Fig-Hi displays the spectrum F{uj) of the normalized ve- 
locity autocorrelation R for two different values of the 
Reynolds number, R c — (100, oo). It is interesting to 
note that there is essentially no inertial range with a 



frequency decay of u> , i.e. 
dependence [Eq. (TIB]) ] , for R e 



a Lorentzian frequency 
= 100; the spectrum de- 
creases as u>~ 4 in the dissipation subrange, i.e. for tur- 
bulent scales close to the Kolmogorov dissipation scale. 

These findings are supported by laboratory exper- 
iments. Fig. [()] displays the measured velocity power 
spectrum in a rotating shear turbulence experiment 
(Karman swirling flow) with water. The solid curve 
displays the measurements and the dot-dashed curve is 
a fit of expression (I17|) to the data with R^J 2 — T^/T^, 
where Tl ~ Te and Ti ~ t n . There is a well de- 
fined inertial range for frequencies between (27tTl) _1 
and (2irT2)~ 1 with a decay of about (w/27r) -2 , i.e. a 
Lorentzian spectrum. The high-frequency tail, how- 
ever, decays more rapidly. The measurement s are su- 
perfici ally similar to the second-order model by Sawfordl 
(|199lh depicted in Fig. [5] It should, however, be noted 
that the molecular Prandtl number of about 6.8 in the 
water experiment is several magnitudes larger than the 
molecular Prandtl number in stellar convection zones, 
w hich is typically of o rder 1 0~ 6 . 



Georgobiani et al.l (|2004l ) used, as did ISamadi et al 



Stci 



<|2003h . the 3D Large-Eddy-Simulations by 
(|200lh to investigate the time-correlation 
Sl(u>, Tfc; r) in the Sun 
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Fig. 9 Comparison of observed (symbols) and simulated 
(curves) stochastic energy supply rates P for solar p modes. 
Simulation results are shown for different sub-grid models 
(from Jacoutot et al. 2008) 



dependence) of the simulation results changes substan- 
tially with wavenumber and depth, and consequently 
that the turbulent energy spectrum function E(k, v) 
is not separable into wavenumber E{k) and frequency 
Q(iv, Tfcjr). The authors suggest using instead the fol- 
lowing empirical function for the temporal part of the 
vertical velocity power spectrum 



1 



(l + A 2 ^ 2 / 



(18) 



function 

Fig.LTI compares the power 
spectrum (squared fast Fourier transform), |w(;y)| 2 , of 
the vertical component of the turbulent velocity field, 
u = (u,v,w), at fixed horizontal wavenumber k and 
depth, between the simulation results (solid curve) and 
various analytical time-correlation functions [Exponen- 
tial (EF), Gaussian (GF) and Lorentzian (LF) fre- 
quency factors] . The results are shown for a wavenum- 
ber k = 4 Mm" 1 computed 250 km below the solar sur- 
face. Neither of the three analytical functions fit the 
simulation results, particularly in the high-frequency 
tail of the spectrum. Moreover, iGeorgobiani et al 



where the two coefficients A and (3 are determined from 
fitting Eq. (TTH)) to the 3D simulation results. Note that 
for /3 = 1 a Lorentzian frequency factor is recovered. 
Results from solar simulations are given in Fig. [5] for 
& Nordlun^ lrl0US hor izontal wavenumbers and d epths. In agree- 

(|2005l ) and iBelkacem et al 



( 2004h reported that the functional form (frequency 



merit with iChaplin et al 

(|2009h the coefficient A oc A increases with depth, i.e. 
the width of the frequency factor decreases with depth 
(see also Fig.^J), particularly at smaller scales (i.e. at 
large wavenumbers). Except near the surface the values 
of the exponent [3 are rather larger than unity. 

4.3 Effect of sub-grid models on simulation results 

Numerical simulations of stellar convection resolve only 
the largest scales of the turbulent motion. The smaller 
scales still need to be approximated by a so-called sub- 
grid model. Such simulations are called Large-Eddy- 
Simulations. The total number of scales that are nu- 
merically resolved, i.e. the ratio of the largest scale Zl 



8 



to the smaller scale l v , is related to the Reynolds num- 

R 3/4 



ber by (e.g. iTennekes fc Lumlevl 119721 ) Il/I v 
and consequently the total number N of grid points 
in the simulation is N = (7l/^) 3 ~ R^ 4 ■ A solar 
Reynolds number of R e ~ 10 12 would require a total 
meshpoint number of N ~ 10 27 . With today's super 
computers the maximum achievable number of mesh- 
points is N max ~ 10 12 , and is therefore some 15 mag- 
nitudes too small for what is required to resolve all the 
turbulent scales of solar convection. Consequently a 
sub-grid model is necessary for describing the dynam- 
ics of the numerically unresolved smaller scales of the 
turbulent cascade. Various models are available. The 
most commonly used models are hyperviscosity mod- 
els and the Smagorinsky model. All sub-grid models 
assume that turbulent transport is a diffusive process. 
Hyperviscosity models, for example, use higher deriva- 
tives for the diffusion operator in the momentum equa- 
tion, thereby extending the inertial range, which also 
leads to a better representation of the dynamics of the 
larger scales. An o verview of sub -grid models was re- 
cently presented by iMieschl (|2005f ). 
An obvious question to ask is what are the effects of 
using different sub-grid models on the mode properties 
in simulations of solar-type sta rs? This quest i on wa s 
recently addressed, in part, bv Uacoutot et alJ (|2008h . 
who compared simulated solar energy supply rates P 
using various sub-grid models. Their results are sum- 
marized in Fig.0 Agreement with observations is gen- 
erally satisfactory except for the classical Smagorinsky 
model. 

It is also important to note that the Prandtl number 
in 3D simulations is currently about 0.01 — 0.25 (e.g 



Miesch et al. 2008). It is therefore substantially larger 



than the Prandtl number in the Sun and in solar-type 
stars. Although numerical simulations provide an im- 
portant input to our understanding of stellar convection 
and are very important for calibrating semi-analytical 
convection models, we must remain aware of the short- 
comings of the currently used 3D numerical simulations. 
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